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Abstract: It is shown that in the noncommutative version of QED (NCQED) Gribov 
copies induced by the noncommutativity of space-time appear in the Landau gauge. This is 
a genuine effect of noncommutative geometry which disappears when the noncommutative 
parameter vanishes. 
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1 Introduction 

Space-time noncommutativity emerges at Plank scale in gedanken experiments when one 
tries to conjugate Quantum Mechanics and General Relativity [1], [2], The result is the 
failure of the classical description of space-time as a pseudo-Riemannian manifold and 
the appearance of uncertainty principles which are compatible with non-commuting co¬ 
ordinates. In different contexts, both string theory [3] and loop quantum gravity [4], [5], 
predict the appearance of space-time noncommutativity at a fundamental level. It is there¬ 
fore natural to investigate how gauge theories have to be modified in order to accomplish 
with the quantum structure of space-time in extreme energy regimes. 

In gauge theory one of the fundamental problems to be solved, in order to be able to 
perform perturbative computations of physical quantities, is the overcounting of the de¬ 
grees of freedom related to gauge invariance (see, for instance, the detailed analysis in [7]). 
The Faddeev-Popov gauge fixing procedure allows for perturbative computations around 
the trivial vacuum A^ = 0. On the other hand, the existence of a proper gauge transfor¬ 
mation preserving the gauge-fixing, would spoil the whole quantization procedure. In [8] 
Gribov showed that in non-Abelian gauge theories (on flat topologically trivial space-times) 
a proper gauge fixing is not possible. Moreover, Singer [9] showed that if Gribov ambigu¬ 
ities occur in the Coulomb gauge, they occur in all the gauge fixing conditions 1 involving 

1 Other gauge fixings are possible such as the axial gauge, the temporal gauge, etc., nevertheless these 
choices have their own problems (see, for instance, [7]). 


- 1 - 





derivatives of the gauge field. The problem is generally addressed in the Landau gauge 
and in this paper we will stick to the latter, although the Gribov-Zwanziger modification 
of the gauge action to cure the problem of Gribov copies has been recently extended from 
the Landau gauge to general R g gauges [10]. 

In the path integral formalism, a Gribov copy close to the identity of the gauge group 
corresponds to a smooth zero mode of the Faddeev-Popov (FP) operator. In order to define 
the path integral in the presence of Gribov copies close to the identity, the most successful 
method is to restrict the path-integral to the neighborhood of A^ = 0 in the functional 
space of transverse gauge potentials, where the FP operator is positive (see, in particular, 
[ 8 ] [11] [12] [13] [14]). When the space-time metric is flat, this approach coincides with the 
usual perturbation theory and, at the same time, it takes into account the infrared effects 
related to the (partial) elimination of the Gribov copies [11] [15] [16]. If one computes the 
propagator corresponding to such a restriction, one finds the famous Gribov form factor 
for the propagator 


G 


G-Z 


(P) 


V 


( 1 . 1 ) 


p A + 7 4 ’ 

where the dimensional constant 7 is related with the size of the Gribov horizon. Although 
at high momenta such a propagator recovers the usual one ~ 1 /p 2 , the infrared behavior 
is drastically different. When one takes into account the presence of suitable condensates 
[17] [18] [19] the agreement with lattice data is excellent [20] [21]. This approach allowed 
to solve (see [22]) the well known sign problem of the Casimir energy and force in the 
MIT-bag model. Thus, in a sense, the Gribov problem is not just a problem since, as 
the whole (refined) Gribov-Zwanziger approach shows, it also suggests in a natural way a 
solution which allows to go far beyond perturbation theory 2 in a very successful way. 

On the other hand, on a space-time with curved metric and/or non-trivial topology the 
situation can be much more complicated since also Abelian gauge theories can have smooth 
zero modes of the Faddeev-Popov operator [23], the maximally (super)symmetric vacuum 
can be outside the Gribov region [24] and the modular region could shrink to zero [25] [26]. 
For these reasons, it is natural to wonder whether the presence of noncommutativity can 
induce Gribov copies even in U(l) gauge theories, which, because of noncommutativity of 
the product, develop self-interaction terms, thus behaving as non-Abelian gauge theories, 
with 


F^ v = d^A u - d v An - i[An, A„]*. (1.2) 

Noncommutativity thus modifies the covariant derivative (see Eq. (2.7) below) in the same 
way as non-Abelian gauge symmetry and it is then natural to expect a non-trivial equation 
for copies. This issue is extremely important since, as we already noticed, noncommutative 
geometry at a fundamental level has been shown to manifest in many different approaches 
to quantum gravity. Consequently, it is mandatory to investigate whether noncommu¬ 
tative geometry induces novel features in QFT which prevent from using the standard 
perturbative techniques. 


2 In particular, both in the case of the glueballs mass spectrum and in the case of (the solution of) the 
sign problem in the Casimir energy and force in the MIT bag model, perturbation theory is obviously not 
sufficient to get the correct answers as non-perturbative physics is needed. 
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The paper is organized as follows. In section 2 we establish the equation for infinites¬ 
imal Gribov copies for noncommutative QED. In section 3 we investigate the existence of 
exact solutions for particularly simple gauge potentials and we actually show that we may 
have an infinite number of genuine noncommutative solutions. In section 4 we discuss our 
results and draw some conclusions. 

2 Equation for Gribov copies 

Let us fix the notations: for each two functions f(x) and g(x) the noncommutative Moyal 
star product / * g is defined as follows 

( f*9)(x ) = f(x)expf^d pff d p d tT \ g{x), (2.1) 

where the indices 3 p, a = 1,.., d and d is the dimension of the space-time, which we assume 
Euclidean. The antisymmetric matrix 9, has the following nonzero components 

01,2 = -02,1 = 01, 03,4 = —04,3 = 02, , 0d-l,d = ~0d,d-1 = 0<i/2, (2.2) 

where 0j are real deformation parameters, in principle all different from each other, charac¬ 
terizing noncommutativity. A rescaling could be performed in order to make all parameters 
equal but we keep them different because it wouldn’t simplify the calculations in the mul¬ 
tidimensional case analyzed in section 3.3. When 0,; —> 0, the star product * goes to the 
standard commutative point-wise product of / and g. 

2.1 Gauge transformation 

Under the 17(1) gauge transformation in NCQED the gauge field A transforms as follows 
A —> A^[a] = U -k Ap-kU^ + iU -k dpU\ U = exp* (ia ), (2-3) 

where the star exponent of an arbitrary function / is by definition 

OO . 

exp*(/) = (2-4) 

n^o n ' 7-- 

n times 

and a is some function of x considered as a parameter of the transformation. It is worth 
noting that in the commutative limit 0 —>• 0, the gauge transformation Eq. (2.3) reduces to 
the standard Abelian gauge transformation 

A —> A^[cc] = Ap + dpQ. + 0(9). (2-5) 

For us of crucial importance will be the infinitesimal form of the gauge transformation 
Eq. (2.3): 

A —> A^[a] = Ap + Dpa + O(a), (2-6) 

where the covariant derivative Dp appears due to non commutativity and is given by 

Dpf = dpf + i (f -k Ap — Ap * /) (2.7) 

for an arbitrary function f(x). 

3 Although it is not necessary, let us assume, that our space time is even dimensional, since we are 
interested in d = 2 and d = 4. 
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2.2 Zero mode equation in general 

Let us fix the gauge to be the Landau gauge, 

= 0 . ( 2 . 8 ) 

In commutative QED this gauge fixing condition fixes the gauge completely, indeed, under 
suitable regularity conditions at the boundary, if A p satisfies Eq. (2.8), the transformed 
field under Eq. (2.5) automatically does not 

d^A'^a] 0, « / 0. (2.9) 

In other words the equation 

d p A^[a\ = 0 (2.10) 

called also the “zero mode equation” has only the trivial solution a = 0. However, in 
a more general setting, for example non-Abelian gauge theories, the zero mode equation 
Eq. (2.10) may have nontrivial solutions a / 0 called zero modes. 

The lack or the presence of zero modes means correspondingly the possibility or impos¬ 
sibility to foliate the functional space of all possible gauge fields into the orbits of the gauge 
group in a way, that each gauge orbit intersects the gauge fixing hypersurface Eq. (2.8) just 
once (in order to be able to integrate over the representatives of each equivalence class). 
The latter situation leads to an overcounting of degrees of freedom when one performs the 
functional integral over A, creating the Gribov problem in the infrared regime. 

The main goal of the present research is to figure out whether the zero mode equation 
Eq. (2.10) can exhibit nontrivial solutions in the case of NCQED, i.e. when A'[a\ is given 
by Eq. (2.3). 


2.3 Infinitesimal Gribov copies 

The infinitesimal zero mode equation which corresponds to the infinitesimal gauge transfor¬ 
mations Eq. (2.6) is of special interest since it has direct relation with the Faddeev-Popov 
ghost action and with the Gribov-Zwanziger term. 

Substituting Eq. (2.6) in the general formula Eq. (2.10) we obtain 


d p D p a = 0. 


( 2 . 11 ) 


Let us understand the structure of this equation from the mathematical point of view. 
Substituting the expression of the covariant derivative Eq. (2.7) and the star product 
Eq. (2.1) into Eq. (2.11) we arrive at the following zero mode equation written in terms of 
a and its derivatives 


- d 2 a + iA p exp - 9 pu d p d a > (d p a) - i(d p a) exp - 9 pa d p d a > A p = 0. 


( 2 . 12 ) 


nonlocal terms 


The presence of nonlocal terms implies that, differently form QCD, this is not a differential 
equation and its resolution is a very hard task. However, in order to say whether we have 
Gribov copies or not we only need to understand whether it has nontrivial solutions a/0. 
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On performing the Fourier transform of a and A one can rewrite the pseudo-differential 
equation Eq. (2.12) as a homogenous Fredholm equation of the second kind. After some 
simple computations we obtain indeed 


0 = -d' z a + iA exp <J O pa d p d a i (d p a) - i{d p a) exp j ^ 8 pa d p d a }> A p 


= I d d ke lkx l — k 2 a(k ) + 2 i I d d q sin I —^0 pa q p k a ) k p A p (q)a(k — (?)}> (2.13) 


1 


which is equivalent to the following integral equation 


k 2 a(k) + 2 i I d d q sin ( -0 pa q p k a ) k p A p (q)a(k — q) = 0. 


Changing the integration variable q —k — q we finally arrive at 

a{k) = / d d q Q(q,k ) a(q), 


(2.14) 


(2.15) 


which is a homogeneous Fredholm equation of the second kind, with the kernel Q given by 

(2.16) 


. 2 ik^AAk — q) ( 1 

Q(g, fc) =-P-— sm f - 9 P q p k a 


1 It is possible to recast the integral equation Eq. (2.15) in such a way that the corresponding 
integral operator becomes manifestly symmetric. To this purpose we notice that, due to 
the gauge fixing condition, A p (k)k p = 0, one can replace k p by (k p + q p )/2 in Eq. (2.16). 
Upon making the change of variable /3 = \k\ - a, with \k\ = ^Jk p k p and multiplying both 
sides of Eq. (2.15) by k, we arrive at the following equivalent integral equation: 


/ 3(*0 = j dd Q P(q,k) /%), 

i (k p + q p ) A p (k - q) . (l opa , 
P(q, k ) =- m -sm 6 p q p k a 


(2.17) 

(2.18) 


In order to show that the linear integral operator P defined by the kernel Eq. (2.18) is 
formally self-adjoint it is necessary and sufficient to show that the corresponding kernel 
satisfies 

P*(q,k) = P(k,q), (2.19) 

where means complex conjugation. Now we recall that we are interested in real gauge 
potentials A(x), which is equivalent to impose 


A*(k) = A(-k) (2.20) 

4 Let us notice that, in the Coulomb gauge djA^ = 0 ,j = 1 1 the equation for the copies is 

formally the same as eq. (2.14), provided we replace k 2 with k 2 and k^A^{q) with k-’Aj(q). Obviously the 
solutions will be different. In particular none of the Ansatze we shall make in the rest of the paper can be 
easily adapted to the Coulomb gauge. We shall comment more on this aspect in the discussion session. 
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on the corresponding Fourier transform. Performing complex conjugation of Eq. (2.18) and 
using Eq. (2.20) and skew symmetry of 9 we immediately obtain 


/i (k^ + q^) A^(k - q ) 

\k\\q\ 



i ( q ^ + k^A^q - k) 

\q\\k\ 



( 2 . 21 ) 


that is exactly the equality Eq. (2.19). 

In principle self-adjoint operators have an infinite set of eigenfunctions and eigenvalues, 
however since we are in the infinite dimensional situation a lot depends on the properties of 
the kernel Eq. (2.18). If for some particular A(k) = B(k), a complete set of eigenfunctions 
ip n with eigenvalues \ n exists, we obtain 


V’n(fc) 


\ r . 


d d qP(q, k)\ A =B^n{q), n= 1,2... 


( 2 . 22 ) 


The latter implies that we have an infinite set of gauge potentials B/X n , n = 1,2... and 
each of them exhibits zero modes a = i^ n . 


2.4 The Henyey approach 

As it is done in standard QCD, also in this case one can follow the Henyey strategy [27] 
where one fixes the form of the zero modes a{k ) and solves for the gauge potential. 

It is worth emphasizing here that, in the standard commutative case, the Henyey 
strategy to fixing the form of the copies and to solving for the gauge potential gives rise to 
algebraic equations which can be easily solved while, in the present nonconunutative case, 
even following such a strategy leads to a rather non-trivial equation, due to the non-locality 
appearing in noncommutative geometry. Therefore we will only sketch the strategy here 
but we will follow a different approach in the next section. 

As we will show, in the noncommutative Henyey approach, instead of a homogeneous 
Fredholm equation of the second kind we obtain a non-homogeneous Fredholm equation of 
the first kind, with Hermitian integral kernel. 

Since A^ has d components constrained by the gauge fixing condition, but the integral 
equation is an equation for one unknown function, we are free to choose some particular 
parametrization of A. Let us make the following Ansatz: 

A^q) = ia(q)G fl (q) withG^q) = -O^q" (2.23) 

with 9 the inverse matrix of the matrix 9. This obviously satisfies the gauge condition 
q^A^{q) = 0. The presence of the imaginary unit in the first line of Eq. (2.23) will become 
clear soon. 

Substituting the Ansatz Eq. (2.23) in Eq. (2.14) we get the following (non homoge¬ 
neous) Fredholm equation of the first kind for the unknown function a(q): 

j d d q R(q, k)a(q) = f(k), (2.24) 
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where 


f(k) = k 2 a(k), and 

R{q,k) = 2 sin Q O^kpq^j k^G p a(k - q). (2.25) 

Note that 

k^Gp(q) = -We^q" (2.26) 

is real and skew symmetric with respect to the exchange k •(-)• q: 

k^Gp(q) = -q^Gp{k) (2.27) 

therefore the combination 

sin Q O pa k p q^j k p G p 

is real and symmetric with respect to the mentioned exchange. Since we are interested in 
real a(.T), the corresponding Fourier transform a(k) satisfies ( a(k ))* = a(—k). Summariz¬ 
ing all observations listed above, we conclude that the kernel R(k, c/)defined by Eq. (2.25) 
satisfies 

(R(k,q))* = R(q,k), (2.28) 

therefore the corresponding linear integral operator R is self-adjoint. Exactly for this reason 
the imaginary unit in the Ansatz Eq. (2.23) is needed. 

In conclusion, in order to solve Eq. (2.24) for the potential it is sufficient to show that 
for some particular choice of a the inverse integral operator R exists, so that a = i? _1 f 
will give us the gauge potential, for which a is a zero mode. 

3 Some exact solutions. 

The question of the existence of Gribov copies is equivalent to the question of the existence 
of eigenvectors of the self-adjoint operator defined by Eq. (2.18). No matter how we choose 
the gauge potential A p , this operator does not belong to the Hilbert-Schmidt class, since 

J dqdk\P(q,k)| 2 = oo, (3-1) 

so one can not say a priori whether there exists at least one gauge potential exhibiting zero 
modes! Neither one can say the opposite. The aim of this section is to demonstrate that 
gauge potentials that give solutions of Eq. (2.18) do exist. 

For the scope of the present section we will not use the property of Hermiticity of 
P, therefore it will be more convenient for the forthcoming computations to resort to the 
original form of the zero mode equation Eq. (2.14), which we rewrite as 

k 2 a(k) + 2ik p J d d q sin (^9 ap q p k a ' S j A p {k — q) a(q) = 0. (3.2) 

We notice that if we consider gauge potentials A p which are proportional to derivatives of 
5{k) , Eq. (3.2) becomes a differential equation for a(k). 
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3.1 The simplest situation 

First we try the following Ansatz 

A tl = (3.3) 

with Q some constant to be fixed. The Fourier transform reads 

A^k) = iQ0^d u 5(k). (3.4) 

This potential obviously satisfies the gauge fixing condition in coordinate space, 

while in momentum space we deal with a distribution, therefore we have to specify in which 
sense the equality 

WA^k) = 0 (3.5) 

holds. Let us fix the space of probe functions a(k) to be the Schwartz space of infinitely 
smooth functions decaying at infinity faster than any arbitrary power. For an arbitrary 
Schwartzian function a(k) one must have 

J d d ka(k)k^A fl = 0 (3.6) 

which is satisfied by Eq. (3.4) it being for arbitrary a(k) 

f d d ka{k)kHQ~9^d v 8{k) = -iQQ^v [d v (k»a{k))\ | fc=Q 

= -iQ(9^a(0) + 6 flu (d"f)\ k=0 • W | fc=0 ) = 0. (3.7) 

o " o " 

The reason to search for zero modes a(x) belonging to the Schwarz space is twofold. On one 
side we observe that, in the commutative case, the zero mode equation is a Laplace equation 
which, unless one specifies boundary conditions, may have nontrivial solutions. Indeed 
each linear function solves it. On the other side the Green function of the corresponding 
Laplacian gives a singular solution, that decreases at infinity for d > 2. In order to get 
rid of these irrelevant solutions (in the commutative case there is no Gribov problem!) we 
impose regularity of a at each finite point and vanishing at infinity, which are both satisfied 
by Schwarz functions. In this class of functions the commutative zero mode equation has 
just the trivial solution a = 0. 

Substituting the ansatz (3.4) in the equation (3.2) and using 

- ‘2QkJ‘6 lw J d d q sin (^0 ap q p k o ^j a(q) q d v 5(k - q) = Qk 2 a(k), (3.8) 

we arrive at the following algebraic equation 

(1 + Q)k 2 a(k) = 0, (3.9) 

which exhibits nontrivial solutions. Indeed if (and only if) 

Q = —1, (3.10) 
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for arbitrary even space-time dimension, any arbitrary function a(k ) is a solution! Unfortu¬ 
nately, although we found nontrivial solutions of Eq. (3.2), this particular gauge potential 
has a peculiar feature. One may show [28] that it is invariant under gauge transformations 
(2.3) and therefore we do not have Gribov copies. 

Nevertheless this potential is of interest. First of all, the existence of such a gauge 
invariant connection is a purely noncommutative feature [28] (also see [29] where such a 
connection has been used to study NCQED as a nonlocal matrix model) and does not exist 
in the commutative limit. Second, its smooth approximations may be used in principle to 
search solutions of the integral equation Eq. (3.2). 

3.2 The next to the simplest situation 

To simplify the presentation let us consider the two dimensional case. Here we have only 
one noncommutative parameter, 9 12 = — #21 = 9. The next to the simplest gauge potential 
leading to a viable differential equation is the following one : 5 

A^(x) oc 9^ v x v x 2 , (3-11) 

which, being in two dimensions, can be further simplified to the form 

A^(x) = Qe^x u x 2 , (3.12) 

with Q some constant to be determined and £^ the Levi-Civita tensor in two dimensions. 
The corresponding Fourier transform reads 

A li {k) = iQe lxv Ud v 5{k). (3.13) 

It is worth emphasizing here that the gauge potential in Eq. (3.13) can be approximated as 
closely as one wants replacing the 5-function with a Gaussian (obviously, if one would use 
the Gaussian from the very beginning the Gribov copies equation would not be solvable 
anymore). Hence, the present example is not only interesting in itself, since it also shows 
that there is a whole family of smooth gauge potentials which are arbitrarily close to having 
smooth normalizable Gribov copies. 

In what follows we will refer to Q as the amplitude of the potential. In spatial coordi¬ 
nates it obviously satisfies the Landau gauge fixing condition and for consistency we check 
whether it satisfies the gauge fixing condition (3.5) in the above mentioned “distributive” 
sense. For an arbitrary probe function a(k ) one obtains 

J ^ka^k^Af, = iQe^ j d d kak )X Ud v 8(k) = - iQe^ [□ d v ( ak :")] | fe=0 

= (Ud u ot) | fe=0 k fl \ k=0 +2(e flu d^d u a)\ k=Q + (Da) \ k=0 e tu/ 6 lu ') = 0. 

' 0 " 0 0 

5 In principle one may also consider quadratic potentials, however for technical reasons we prefer to 
deal with rotationally invariant potentials, therefore the next to the simplest potential which we consider 
is cubic. Indeed, in d, = 2 rotationally invariant quadratic in x gauge fields do not exist. 
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Let us now substitute the potential Eq. (3.13) in the integral equation Eq. (3.2) in order 
to derive a partial differential equation for the zero modes a(k). We obtain 


J d a q ( 9 D q d v 5(q - k )) sin ( ^9 ap q p k a ) a(q) 

= -Qk p e pv Un q d v 

Q9 


sin ( -9 ap q p k a ) a(q) 


q=k 


0 2 k 4 a — 4 fc 2 Dd — 8 s pu e r,x k lJ ,k rt d u d\d^ (3-14) 

hence the zero modes d(A:)have to satisfy the partial differential equation given below: 

4fc 2 D — 8 e pv £ r,x k p k ri d u d\ — a(k) = 0. (3.15) 

We notice that, since in two dimensions e pu is a universal tensor, this equation is rota- 
tionally invariant, therefore it makes sense to rewrite it in polar coordinates (r, <p) given 
by 

{ k\ = r cos 4> 
k 2 = r sin (j) 


(3.16) 


One may easily see that in polar coordinates Eq. (3.15) reads 

r 2 a rr + 3rb r + ~ “j“ r4 “ + = 0, 


(3.17) 


therefore it exhibits separation of variables. Let us look for a solution in the following form 

<*(<!>, r) = 4>(0)/(r), (3.18) 

where the functions and / satisfy the following ordinary differential equations 


-<$> H = A4>, 

r 2 f rr + 3 rf r + 3A + ^r 2 - ^r 4 ^) / = 0. 


(3.19) 


The former is just the equation of the simple harmonic motion, while the latter is the 
confluent hypergeometric equation whose properties are very well studied, see e.g. [31] for 
a review. Let us specify the boundary conditions as follows: 


4>(0) = 4>(27t), 

|/(0)| < oo, /(r) -» 0, whenr -> oo. 


(3.20) 


Below we will see that each function satisfying the boundary conditions given by Eq. (3.20) 
belongs to the Schwarz space. We also notice that the deformation parameter 6 enters the 
equation Eq. (3.19) in the combinations Q9 and 9 2 , where Q is arbitrary. Therefore without 
loss of generality one may consider only 9 > 0, since the opposite sign of 9 corresponds to 
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the opposite sign of the arbitrary amplitude Q. From the angular boundary conditions we 
see that 

A n = n 2 , n = 0, ±1, ±2,... (3.21) 

so that the general solution for the angular equation Eq. (3.19) is of the form 

<3?(0) = ci cos (ncj>) + C 2 sin (ncj)), (3.22) 


where ci and C 2 are arbitrary constants. The general solution of the radial equation for 
A = An is given by 


/(r ) = r vra-i eX p 

y y _ y 

with a=- + -v'3n2 + l-^, 


Or* 


6 r 2 


-— ) [ciM +c 2 U {a,c,— 

c = 1 + y/ 3n 2 + 1 


(3.23) 


U and M are Rummer functions and ci and C 2 are arbitrary constants. 

We notice however that the boundary condition for the radial dependence can be 
satisfied if and only if the number a defined by Eq. (3.23) is a non positive integer 

a = — m, m = 0,1,2,... (3.24) 

In that case Rummer functions reduce to Laguerre polynomials. 

Thus, the solution regular both at zero and infinity (see [31]) is 

/( r ) — Cr v w+t-i ex p (-d^, (3.25) 

where C is an arbitrary constant, and L“(,z) stands for generalized Laguerre polynomial. 
This solution exists when the amplitude Q takes one of the discrete values 

1 


Qnm — 


n = 0, ±1, ±2,..., m = 0,l,2,... 


(3.26) 


0 2 (y/3n 2 + 1 + 2m + 1)' 

The general form of the zero modes, when the amplitude Q belongs to the discrete set 
defined above, is 

a n m(r, 0) = (Ci cos {ncj)) + C 2 sin (rid))) C /3n2+1 ~ 1 exp 
where Ci, C 2 are real if n is even 

and Ci,C 2 are purely imaginary if n is odd. (3.27) 


The restriction on arbitrary constants comes from the requirement to have real a(x). In¬ 
deed, in order to satisfy this requirement, the corresponding Fourier transform a must 
satisfy ( a(k ))* = a(—k). From another side the reflection k — > —k is equivalent to the shift 
4> —> <p + vr in polar coordinates. One may easily check that the radial dependence is real, 
and the linear combination of sine and cosine appearing in (3.27) satisfies 

(Ci cos ( n((f> + vr)) + C 2 sin (n(0 + 7r))) = (Ci cos (ncj)) + C 2 sin (ncj)))* (3.28) 

if and only if the restriction described in (3.27) is imposed. 
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l 

Pi 

Pi 

1 

0 

1 

2 

1 

2 

3 

4 

7 

4 

15 

26 

5 

56 

97 

6 

209 

362 





Table 1. The table shows the first six numbers from the A001353 sequence. 

3.2.1 How many zero modes do we have? 

In this subsection we discuss how many linearly independent solutions for a given ampli¬ 
tude (3.26) we have, connecting our results with number theory. From the solution (3.27) 
we can guarantee at least two. It is remarkable that for some special subset of the am¬ 
plitudes Q n m we can have more than two, or more precisely the following statement holds: 

Theorem. For arbitrarily large number N there exists such an amplitude Qnm, that the 
number of linearly independent zero modes corresponding to this potential is greater than 
N. 

Before we go ahead to prove the existence of “many” solutions, let us remind some facts 
from number theory. The key property allowing such an existence is due to the fact that 
there exist infinitely many natural numbers pi, l = 1, 2,... such that 

3 pf + 1 is a perfect square, (3.29) 

i.e. there exists such a natural number Pi, that 

llpf + \ = P?, 1 = 1 , 2 ,... (3.30) 

This sequence, called “A001353”, is well known and studied (see [30] and refs therein). We 
have 

Pi = 0, P2 = 1, Pi = tyi-i -Pi-2, 1 = 3,4,... (3.31) 

and the corresponding Pi are given by 

P l = 2p l -pi_ 1 . (3.32) 


We also notice that all P 2 k , fe = 1,2,... are even natural integers, therefore 


Jk = ~y, k = 1,2,... (3.33) 

are natural integers. The first six numbers from the sequence A001353 are presented in 
Table 1 . 
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k 

m 

V3n 2 + 1 

n 

1 

Jn - 1 

2 

i 

2 

Jn — 13 

26 

15 

3 

J N ~ 181 

362 

209 





k 

Jn — Jk 

Pk 

P2k 





N 

0 

Pn 

P2N 


Table 2. This table summarizes the construction of N different solutions for a given amplitude, as 
described in section 3.2.1. 


We are now ready to prove the above stated theorem. For a given (arbitrary large) N, 
we choose the amplitude Q = Qn where 


Qn = 


1 


e 2 (2Jjv +1)' 

This implies that the parameter a appearing in Eq. (3.23) 


(3.34) 


_ 1 ( V3n 2 + 1 1 

" “ 2 H 2 WQ 


— -Jn + 


V 3n 2 + 1 


(3.35) 


is a nonpositive integer if and only if ^ 3 ^ 2+1 i s a positive integer not greater than Jn- Let 
us now substitute n = ri}~ given by 


n k = P2k j A: = 1,..,W (3.36) 

We obtain 

a = — Jn + Ju, k = l,...,N, (3.37) 

and for each k from 1 to N we have two linearly independent solutions (with sine and 
cosine), given by (3.27), where the Laguerre polynomial is labelled by m = —Jn + Jk-, 
and the parameter \/3n 2 + 1 = P}~. Finally we conclude that we have found the amplitude 
Q = Qn for the gauge potential in (3.12), that exhibits 2N linearly independent zero 
modes. QED 

The logic of our construction is summarized in Table 2. 

3.2.2 Examples 

Let us illustrate the results of the previous section with some examples. 


• N=l. The corresponding amplitude Q is given by (see Eq. (3.34)) 

Ql = 7^2- 


(3.38) 
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There are just two linearly independent solutions 


m = 0, n = 1 angular dependence = cosine 
m = 0, n = 1 angular dependence = sine 
In Fig. 1 we plot the corresponding potential for the cosine case 6 . 



k 


i 


Figure 1. Shape of ia, for N=1 and angular dependence of cosine type. The deformation 
parameter 9 is chosen to be equal to one. 

• N=2. The corresponding amplitude is given by (see Eq. (3.34)) 

= 2702 ' 

There are four linearly independent solutions 


m 

= 12, 

n = 1 

angular dependence = 

cosine 

m 

= 12, 

n = 1 

angular dependence = 

sine 

m 

= 0 , 

n = 15 

angular dependence = 

cosine 

m 

= 0 , 

n = 15 

angular dependence = 

sine 


In Fig. 2 we plot the corresponding potential with angular dependence of cosine type. 


6 Since for odd n the Fourier image a(k) is purely imaginary, we multiplied it by i to build a plot. 


(3.39) 
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Figure 2. Shape of ia, for N=2. The first plot represents the case to = 12, while the second one 
corresponds to to = 0 of cosine type. The deformation parameter 9 is chosen to be equal to one. 


3.3 Multidimensional generalization 

Let us consider now the generalization to arbitrary dimensions of the solutions discussed 
above. 

Although the potential given by Eq. (3.11) satisfies the Landau gauge fixing condition 
in arbitrary even dimensions, it does not lead to a simple partial differential equation where 
one can easily separate variables. However, if one considers the space as a direct product 
of d/2 orthogonal planes, the tensor e / “' is invariant under rotation in each plane. This 
observation suggests to define the potential by reproducing the two dimensional structure 
in each plane: 

A\{x) = Qix + x%) 

A 2 {x) = —Qix i(xf + x%) 

< A 3 (x) = Qnx 3 (x§ + x\) (3.40) 

A 4 (x) = -Qnx 4 (x| + x'j) 


where Qi,Qn ,... are in general different constants. 

Taking the Fourier transform of Eq. (3.40) and substituting the result in the integral 
equation (3.2), after carrying out similar computations to the ones that we did in the d = 2 
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case, we arrive at the following partial differential equation 


[QiDj + QiiDii + •••] & — (fcj + kjj + (3.41) 

where we used the following notations 

kj = kf + k %, 
kji = + k'l , 

. (3.42) 

and 

Dj = -0i (d? + df) + 2 k\dl + 2 k 2 2 dl - Ak x k 2 d 2 vl - 6 -j-ki^j 

D n = 9 2 (k 2 rr (<9f + <9f) + 2&f<9f + 2kjd 2 - 4 fc 3 fc 4 9| 4 - 

. (3.43) 

with Oi, i = l,..., d/2 the different noncommutative parameters of each two-dimensional 
plane, as in Eq. (2.2). One may easily see that Eq.(3.41) is a sum of d/2 equations 

QiDja = akj , (3.44) 

QnDna = ak 2 j, (3.45) 

. (3-46) 


where each equation is of the form (3.15). 

If functions ai(ki, k 2 ), ajj(k 3 , Aq),..., are solutions of two dimensional equations Eq. (3.44), 
Eq. (3.45), ..., then their product 

& = &i(ki, k 2 ) ■ au(k 3 , fc 4 ) • ... (3.47) 

solves each of the equations Eq. (3.44), Eq. (3.45), ... as well as their sum Eq. (3.41). 

But we have already constructed solutions for the two dimensional case belonging to the 
Schwarz space 5(M 2 ) hence the product Eq. (3.47) automatically belongs to the Schwarz 
space S(R d ). 

4 Discussion and outlook 

In the present paper we have shown that an infinite number of Gribov copies exists in 
noncommutative QED, and this is a genuine noncommutative effect. As already recalled 
in the introduction and well known in the literature on noncommutative gauge theory, 
NCQED behaves like a non-Abelian gauge theory because of the non-trivial deformation 
of the covariant derivative Eq. (2.7). Therefore the existence of a nontrivial equation 
for copies was expected and not surprising. The main result of the paper is that such 
equation has solutions, which we compute, and they are an infinite number. The intrinsic 
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interest of this observation is that it shows how noncommutative geometry can give rise 
to a global obstruction preventing a proper gauge-fixing already in Abelian gauge theories. 
Consequently, the intriguing possibility which is naturally suggested by the present analysis 
is to extend to the case of NCQED the Gribov-Zwanziger restriction. Indeed, the Gribov- 
Zwanziger restriction would yield to the following modification of the NCQED propagator 


G G ~ Z (p ) 



(4.1) 


with 7 depending on the noncommutative parameter 0. 

However, having proved that the noncommutativity of space-time induces (infinitely 
many) Gribov copies also in NCQED (as we did in the present paper) is not enough to 
justify the Gribov-Zwanziger approach. 

In non-Abelian gauge theory on flat topologically trivial space-times, the Gribov- 
Zwanziger approach is based on the following fundamental results [12]: 

1) The Gribov region is bounded in every direction (in the functional space of trans¬ 
verse gauge potentials). 

2) The Faddeev-Popov determinant changes sign at the Gribov horizon. 

3) Every gauge orbit passes inside the Gribov horizon. 

The importance of the last result lies in the fact that it justifies the restriction of 
functional integration to the Gribov region (since what is left outside the Gribov horizon 
is just a copy of something inside it and so no relevant configuration is lost). 

In order to justify the Gribov-Zwanziger restriction in NCQED we should generalize 
the analysis of [12] to the noncommutative case. This is a highly non-trivial technical 
task since many of the arguments used in these references to prove the properties 1 ), 
2) and especially 3) make heavy use of the theory of local elliptic PDEs while, in the 
noncommutative case, the Gribov copy equation becomes non-local. We hope to come 
back in a future publication on these important issues. 

As a final remark, it is natural to wonder whether the present results can be extended 
to the Coulomb gauge as well. In the footnote 4 we have described the equation of the copies 
in such a case. Formally, some of the copies constructed here can also be used to construct 
copies in the Coulomb gauge (for instance, a copy in the Landau gauge in a two dimensional 
non-commutative space-time can be trivially promoted to a copy in the Coulomb gauge in 
3+1 dimensions). However, the physical interpretation is rather subtle. Indeed, one of the 
fundamental properties of the Moyal star product is the ‘democracy’ between all space-time 
coordinates. Obviously, such democracy is not respected by the Coulomb gauge (one well 
known consequence being the existence of residual gauge transformations in this gauge). 
Thus, in the non-conunutative setting, not only one should declare the Euclidean time as 
‘special’ but its non-commutative partner as well and this would make the analysis much 
more complicated than in the Landau case. We hope to come back on this interesting issue 
in a future investigation. 
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